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Abstract 

The lifetime and oscillation frequency of the B° meson has been measured using B° — > D* + l~v 
decays recorded on the Z° peak with the OPAL detector at LEP. The D* + — ► D°7r + decays were 
reconstructed using an inclusive technique and the production flavour of the B° mesons was deter- 
mined using a combination of tags from the rest of the event. The results 

r B « = 1.541 ± 0.028 ± 0.023 ps, 
Am d = 0.497 ± 0.024 ± 0.025 ps" 1 

were obtained, where in each case the first error is statistical and the second systematic. 
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1 Introduction 

The lifetimes of b hadrons depend both on the strength of the b quark coupling to the lighter c and 
u quarks, and on the dynamics of b hadron decay. The spectator model prediction that the lifetimes 
of all heavy hadrons containing the same heavy quark are equal is modified by non-spectator effects 
dependent on the flavour of the light quark(s) in the hadron. In contrast to the charm hadrons, where 
t d + sa 2.5t£)0 01 j non-spectator processes are expected to lead to lifetime differences between the 
B + and B° mesons of at most 10% [||. Measurements of b hadron lifetimes at the level of a few 
percent or better are therefore needed to test these predictions, and probe the non-spectator processes 
contributing to the decays. In addition, precise measurements of the B° lifetime are also needed for 
the determination of the magnitude of the CKM matrix element V c b Pj • 

The most precise measurements of the B + lifetime come from topological vertex reconstruction 
techniques, where the selection of charged secondary vertices allows a clean sample of B + decays to be 
isolated J§,f|. This method is however limited for B° lifetime measurements, due to contamination from 
other neutral b hadrons (B s mesons and b baryons). An alternative technique is to use B° — > T)* + l~9 
decays^, which can be efficiently partially reconstructed by exploiting the low energy release in the 
decay D* + — ► D°7r + . In this case, only the ir + from the D* + decay is identified, and no attempt is 
made to fully reconstruct the D° meson decay. This method has previously been used by DELPHI || 
to measure the B° lifetime, and by DELPHI and OPAL || to measure |Tcb|- 

The same sample of B° — ► Y)* + l~v decays can also be used to measure the B°-B° oscillation 
frequency^ Amj |J. In the neutral b meson system, the weak eigenstates B° and B° differ from 
the mass eigenstates, and transitions between them are possible, arising dominantly in the Standard 
Model from second order weak transition box diagrams involving virtual top quarks. Therefore an 
initial B° meson can oscillate into a B° at time t with a probability given by 

P(B° -» B°) = -(1 - cosAm d i) . 

Measurements of Am^ allow the extraction of the magnitude of the CKM matrix element Vtd, though 
the precision is currently severely limited by theoretical uncertainties |lOf| . 

In this paper, measurements of both t^o and Am^ based on this technique are presented. The 
reconstruction of B° — > D* + £~v decays is described in Section |2[ followed by the determination of 
the proper decay time for each event in Section |3], and the production flavour tagging needed for the 
oscillation measurement in Section |4|. The fit to determine r B o and Am^ is described in Section ||, 
followed by a discussion of systematic uncertainties in Section ||. The results are summarised and 
combined with previous OPAL measurements in Section |?]. 



2 Inclusive reconstruction of B° D* + £-p events 

The OPAL detector is well described elsewhere [11]. The data sample used in this analysis consists of 



about 4 million hadronic Z° decays collected during the period 1991-1995, together with an additional 
400 000 events recorded primarily for detector calibration purposes in 1996-2000. Corresponding 
simulated event samples were generated using JETSET 7.4 [^] as described in Jl3|] . 

Hadronic Z° decays were selected using standard criteria |l]|. To ensure the event was well 
contained within the acceptance of the detector, the thrust axis direction^] was required to satisfy 
| cos#t| < 0.9. Tracks and electromagnetic calorimeter clusters with no associated tracks were then 
combined into jets using a cone algorithm WM, with a cone half-angle of 0.65 rad and a minimum jet 



1 Charge conjugate reactions are implied, and the symbol I refers to either an electron or muon. 
2 The conventions c — 1 and U = 1 are employed throughout. 

3 A right handed coordinate system is used, with positive z along the electron beam direction and x pointing to the 
centre of the LEP ring. The polar and azimuthal angles are denoted by 6 and <j>. 
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energy of 5 GeV. The transverse momentum pt of each track was defined relative to the axis of the jet 
containing it, where the jet axis was calculated including the momentum of the track. 

The reconstruction of B° — > D* + l~is events was performed by combining high p and pt lepton 
(electron or muon) candidates with oppositely charged pions from the D* + — > D°7r + decay. The 
selection is similar to that used in || , but with some changes to produce an unbiased B° decay proper 
time measurement and to increase the efficiency at the expense of higher combinatorial background. 



Electrons were identified and photon conversions rejected using neural network algorithms [13|, and 
muons were identified as in |||. Both elect rons and muons were required to have momenta p > 2 GeV, 
transverse momenta with respect to the jet axis pt > 0.7 GeV, and to lie in the polar angle region 
|cos0| < 0.9. 

The selection of pions from D* + decays relies on the small mass difference of only 145 MeV |jj be- 
tween the D* + and D°, which means the pions have very little transverse momentum with respect to 
the D° direction. In each jet containing a lepton candidate, the D° momentum vector p D o and energy 
were reconstructed using an inclusive technique, selecting a group of tracks and calorimeter clus- 
ters compatible with a D° decay using kinematic, impact parameter and invariant mass information. 
This procedure is described fully in || and gives angular resolutions on the D° direction of about 
45 mrad in both 6 and 0. 

Each track in the jet (other than the lepton) was then considered in turn as a slow pion candidate, 
provided it satisfied 0.5 GeV < p < 2.5 GeV and had a transverse momentum with respect to the D° 
direction of less than 0.3 GeV. If the pion under consideration had been included in the reconstructed 
D°, it was removed and the D° momentum and energy recalculated. The final selection was made 
using the reconstructed mass difference^ AM between the D* + and D° mesons, calculated as 



AM = y/£%. - |p D *| 2 -m B o , 

where the D* + energy is given by E-q* = E D o + E w and momentum by pd* = Pd° + Pn- 

The position of the B° candidate decay vertex was reconstructed from the intersection point of 
the lepton and slow pion tracks in the x-y plane. The two-dimensional flight distance of the B° was 
then calculated as the length of the vector between the e + e~ interaction point ('beamspot') and the 
B° decay vertex, constrained to lie along the x-y projection of the jet direction. This was converted to 
a three dimensional decay distance L using the jet polar angle. Using just the lepton and pion tracks, 
together with the interaction point position and uncertainty determined with a fit to many consecutive 
events , results in a decay length estimate which is bias free and whose resolution does not depend 
strongly on the decay length itself. Although the decay length resolution could be improved by adding 
more tracks to the B° decay vertex, this would introduce significant bias at small decay lengths, and 
is not necessary as the resolution is already adequate for the measurement of r B o and Awid. 

The reconstructed decay length L was signed positive if the B° decay vertex was displaced from 
the beamspot in the direction of the jet momentum, and negative otherwise. The decay length error 
cfl was calculated from the track parameter and beamspot position error matrices. The decay length 
and error were required to satisfy —0.5 cm < L < 2 cm, L/ai > —3 and ul < 0.2 cm. 

The resulting distributions of AM for opposite and same sign lepton-pion combinations are shown 
in Figure |]( a) and (b). The predictions of the Monte Carlo simulation are also shown, broken down 
into contributions from signal B° — > D* + £~P events, 'resonant' background containing real leptons 
combined with slow pions from D* + decays, and combinatorial background, made up of events with 
fake slow pions, fake leptons or both. 

In Monte Carlo simulation, about 36% of opposite sign events with AM < 0.17 GeV are signal 
go events, 15 % are resonant background and 49 % are combinatorial background. The 

resonant background is made up mainly of B — > D* + ir~£~i>, B° — > Y)* + -K®l~i> and B s — > D* + K°^~P 
decays. These are expected to be dominated by b semileptonic decays involving orbitally excited 



4 The D* + -D° mass difference AM was denoted by Am in 
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Figure 1: Reconstructed AM distributions for selected (a) opposite sign and (b) same sign lepton-pion 
combinations. The data are shown by the points with error bars, and the Monte Carlo simulation con- 
tributions from signal B° — > ~D* + t~i> decays, other resonant D* + decays and combinatorial background 
are shown by the open, single and cross hatched histograms respectively. 
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charm mesons (generically referred to as D**), e.g. B — > D**°£~v followed by D**° — > D*+7r~. These 
decays will be denoted collectively by B — > D* + 1i£~za Small contributions are also expected from 
b — > D* + r^X decays (via any b hadron) with the r decaying leptonically, and b — > D* + D~X with 
the D~ decaying semileptonically (each about 1 % of opposite sign events). For same sign events 
with AM < 0.17 GeV, there is a small resonant contribution of about 6% from events with a real 
D* + — > D°7r + where the D° decays semileptonically, and the rest is combinatorial background. 



3 Proper time reconstruction 

The proper decay time t of each B° candidate was calculated from its reconstructed decay length 
L and energy E-qo. The energy was calculated using a technique similar to that described in [|17|], 
exploiting overall energy and momentum conservation in the event to account for the missing energy 
of the unreconstructed neutrino. The event was treated as a two-body decay of a Z° into a B° of mass 
5.279 GeV M and another object making up the rest of the event. The B° energy was calculated as 



B° 



^cm + m B° ~~ ^rec 



2-E/cm 



where E cm is the centre-of-mass energy of the event and M rec the invariant mass of the object recoiling 
against the B°. The latter was calculated from all tracks and calorimeter clusters in the event, 
excluding the lepton and those associated to the reconstructed D* + . A correction for double counting 



of charged particles in the tracking detectors and calorimeters was applied [18|, and the recoil mass 
was first scaled by {E/E V - 1S ) where -B v i s is the total event visible energy and E = 87 GeV is the typical 
visible energy in events with only one neutrino. This procedure improves the resolution in events 



where a second neutrino is present [18]. The resulting energy estimate is unbiased and has an RMS 
resolution of 3.8 GeV in Monte Carlo B° — > D* + £~i> events, as shown in Figure ||(a). 

The proper decay time t was then calculated from the candidate decay length and energy as 

m B oL 



El 



B o m B o 

The resulting proper time resolution depends on the true proper time t 1 , but has no significant bias, 
as shown in Figure |2|(b-f). The resolution degrades with increasing proper time due to the larger 
influence of the energy resolution at large decay lengths. The resolution was parameterised as a 
function i? D .+ (£,£') giving the expected distribution of reconstructed t for each true value t'. The 
resolution function was implemented as the sum of three Gaussian distributions, whose widths vary 
linearly with t'. The resolution function is also shown as the solid line in Figure ^(b-f), and gives a 
reasonable description of the Monte Carlo resolution, adequate for measuring the lifetime and relatively 
slow B° meson oscillations. 

Similar resolution functions Rj)*+ n - (t, t'), R^*+^o(t,t') and i2 D *+ K o(i, i') were generated for the 
three main D** background contributions B — > D* + ir~l~v, B° — ► D* + ir°£~v and B s — > D* + K°£~v. 
In these events, the b hadron energy reconstruction has a bias of about —1.5 GeV as some of the b 
hadron decay products are incorrectly included in the recoil mass. This bias was corrected in the 
corresponding resolution functions. 



2 



4 Production flavour tagging 

The B° oscillation frequency measurement requires that the production and decay flavour (B or B°) 
of each meson be determined, in addition to the decay proper time t. The decay flavour can be 
determined from the sign of the lepton in the B° — > ~D* + £~v decay, but the production flavour must 
be determined from other information in the event. The production of quark anti-quark pairs in 
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Figure 2: (a) Reconstructed b hadron energy resolution and (b-f) reconstructed proper time resolution 
in various ranges of true proper time t', for Monte Carlo B° — > T>* + £~u events. The Monte Carlo is 
shown by the points with error bars, and the resolution function is shown by the solid line in (b-f). 
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Z° — » bb decays allows the production flavour of the B° to be inferred from that of the other b hadron 
in the event. This information is diluted due to the possible mixing of the second b hadron. The 
other particles in the B° jet (produced in the fragmentation of the b quark) also carry some useful 
information. 

The tagging information was extracted using the methods described in 20 1. The event was 



divided into two hemispheres by the plane perpendicular to the thrust axis and containing the e + e~ 
interaction point. In the hemisphere opposite to that containing the B° jet, up to three pieces of 
information were used to tag the b hadron flavour: 

• The jet charge Q opp of the highest energy jet in the opposite hemisphere, defined as 

n _ Ei(Pi) K qi m 
Qopp ~ ZM K [ ) 

where p\ is the longitudinal momentum component with respect to the jet axis and qi the charge 
of track i, and the sum is taken over all tracks in the jet. The parameter k was set to 0.5, as 
in |0|. 

• The charge Q v tx of a secondary vertex (if existing), reconstructed in any jet in the opposite 
hemisphere, as in 0|. A well reconstructed charged vertex indicates a B + or B hadron, tagging 
the parent quark as a b or b respectively, whilst a neutral or badly reconstructed vertex gives 
no information on the b quark flavour. 

• The charge of a high momentum lepton found in any jet in the opposite hemisphere, selected as 
in Q . A high momentum lepton is most likely to come from a b hadron decay, again tagging the 
parent b or b quark according to its charge. A neural network algorithm was used to suppress 
fake leptons and those coming from cascade charm decays (b — > c — > £) which have the wrong 
charge correlation. 

These variables were combined using a neural network algorithm into a single tagging variable Qt 
for the opposite hemisphere @,|2(|, quantifying the confidence with which the hemisphere was tagged 
as containing a b or b hadron. Different neural networks were used depending on what combination 
of vertex and /or lepton variables were available to combine with the jet charge Q opp . 

In the hemisphere containing the B° jet, only the jet charge Q sa me can be used to infer the B° 
production flavour. Qsame was calculated using equation ||, but with the parameter k set to zero, so it 
becomes simply the average of the charges of the tracks in the jet. This avoids being sensitive to the 
decay flavour of the B° (and hence whether it has mixed or not), but is still sensitive to the production 
flavour via the information carried by the fragmentation tracks in the jet . The jet charge Q S ame was 
used to generate a second hemisphere tagging variable Qm independent of the opposite hemisphere 
variable Qt- 

The two tagging variables were combined to produce a single tag Q2 for the B° production flavour, 
as in [||. The continuous variable Q2 ranges from —1 to +1, and is defined such that events with 
Q2 = +1 are tagged with complete confidence as containing a produced B°, events with Q2 = — 1 are 
tagged with complete confidence as containing a produced B°, and events with Q2 = are equally 
likely to be B° or B°. The modulus IQ2I satisfies IQ2I = 1 — 277, where rj is the 'mis-tag' probability, 
i.e. the probability to tag the production flavour incorrectly. Finally, the production flavour tag Q2 
and B° decay lepton sign / were combined to produce the mixing tag Q = Q2 ■ I, such that events with 
Q > (Q < 0) are tagged as unmixed (mixed). 

Both the sign and magnitude of Q are used in the fit to determine Amj, giving the events with 
high probability to be correctly tagged more weight. Considering only the sign of Q, 31 % of signal 
events are tagged incorrectly. The event-by-event weighting reduces this to an effective mis-tag of 
28 %, equivalent to a 33 % increase in statistical sensitivity. 
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5 Fit and results 

The values of t^o and Am^ were extracted using an unbinned extended maximum likelihood fit to 
the reconstructed mass difference AM, proper time t and mixing tag Q of each event. Both opposite 
and same sign events with AM < 0.22 GeV and —2 < t < 15 ps were used in the fit, the high AM 
and same sign events serving to constrain the combinatorial background normalisation and shapes in 
the opposite sign low AM region populated by the B° — > Y)* + £~v decays. Using the AM value from 
each event in the fit, rather than just dividing the data into low AM 'signal' and high AM 'sideband' 
mass regions, increases the statistical sensitivity as the signal purity varies considerably within the 
low AM region. 

The likelihood is similar to that used in [|| for the measurement of |V^b|- The logarithm of the 
overall likelihood was given by 

M a M*> 

ln£ = ^ln£? + ^ln4-iV a -iV fe (2) 
i=l j=l 

where the sums of individual event log-likelihoods In £f and In C b are taken over all the observed M a 
opposite sign and M b same sign events in the data sample, and N a and N b are the corresponding 
expected numbers of events. 

The likelihood for each opposite sign event was given in terms of different types or sources of event 

by 

4 

£"(AMj, U, Qi) = E NgM s (AMi) T s (t u Q { ) (3) 

s=l 

where iV° is the number of expected events, M S (AM) the mass difference distribution and T s (t,Q) 
the proper time distribution for source s. For each source, the mass difference distribution M S (AM) 
is normalised to one. The total number of expected events is given by the sum of the individual 
contributions: N a = J2t=i N£. 

There are four opposite sign sources: (1) signal B° — > D* + l~i> events, (2) B — > D* + h£ _ z/ events 
where the D* + is produced via an intermediate D**, (3) other opposite sign background involving a 
genuine lepton and a slow pion from D* + decay and (4) combinatorial background. The sum of sources 
2 and 3 are shown as 'resonant background' in Figure |l[ A similar expression to equation |3| was used 
for C'j, the event likelihood for same sign events. In this case, only sources 3 and 4 contribute. 

The mass difference distributions M S (AM) for sources 1-3 were represented by analytic functions, 
whose parameters were determined using large numbers of simulated events. For the signal (source 
1), both unmixed and mixed events must be considered to determine the proper time distribution 
Ti(t, Q). The probability to find a signal decay with true proper time t' is given by: 

Pi(t',Q) = J-e-V-BO ( P U (Q) C 1 + «* + P m (Q) ( 1 ~ «" W) 

where P u (P m ) is the probability that the event is unmixed (mixed) given the observed mixing tag Q. 
These probabilities are given by P U (Q) = (1 + Q)/2 and P m (Q) = (1 — Q)/2. The signal probability 
is then convolved with the time resolution function i? D *+(i,t') described in Section |3] to give the 
expected reconstructed proper time distribution T\(t, Q) as a function of the assumed B° lifetime t-qo 
and oscillation frequency Amj: 

roo i 

Ti(i,Q)=/ dt'- e-* /r B°(l + QcosAm d O^D*+(t,0- (4) 

Jo 2r B o 

The corresponding proper time distribution T2(t, Q) for source 2 (D**) was calculated from the sum of 
the individual contributions from B + , B° and B s decays via D**. For B° and B s decays, equation was 
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Quantity 


Assumed value 


Reference 


T B+ 


1.65 ±0.03 ps 


M 




1.54 ± 0.07 ps 


1 


Am s 


9.1-50 ps" 1 


p2[. see text 


Rb 


(21.70 ± 0.09) % 




Br(D*+ -> D°7r+) 


(68.3 ± 1.4)% 


PI 


Br(b -> D*+h£i/) 


(0.76 ±0.16)% 


id 


Br(b -> D*+r-pX) 


(0.65 ±0.13)% 


11 


Br(B° D*+D^ )_ ) 


(4.2 ± 1.5) % 


i 


Br(b -> D*+X) 


(17.3 ± 2.0) % 




Br(c -f D*+X) 


(22.2 ±2.0)% 


pi 



Table 1: Input quantities used in the fit for tbo and Amj. 



modified, with the resolution function i? D ,+ (i,t') replaced by Ri)*+ n o (t, t') (B° decays) or i? D *+ K o(t, i') 
(B s decays), and rgo and Amj replaced by tb s and Am s , the B s lifetime and oscillation frequency, for 
B s decays. For B + decays, which are always unmixed, the proper time distribution takes the form 

poo 1 

T 2 , B+ (t, Q) = / dt'- e ~* / r B+ (1 + Q)i2 D . +7r - (t, t') . 

Jo 2r B + 

The values of Tg+ , tb s and Am s and their corresponding uncertainties were taken from |l], ^] and are 
given m Table Q Ara s was varied between the experimental lower limit and 50 ps , well above the 
experimental sensitivity and theoretically favoured region. 

The number of signal events Nf was left as a free parameter in the fit. The number of events in 
source 2, and the relative contributions of B + , B° and B s decays within this source, were calculated 
as described in H, from the measured branching ratio Br(b — ► D* + 7r~^PX) = (0.473 ± 0.095) % pjj 
and the assumption of isospin and SU(3) symmetry. This was combined with the selection efficiency 
estimated from Monte Carlo and the values of and Br(D* + — > D°7r + ) listed in Table || to determine 
the expected number of these events in the data sample. 

The numbers N^ b of events in the small background contributions covered by source 3 (both 
opposite and same sign) were taken from Monte Carlo simulation, with branching ratios adjusted to 
the values given in Table ^ as described in more detail in ||. The proper time distributions T-^(t, Q) 
were described by negative exponentials convolved with Gaussian resolution functions, with parameters 
again determined from Monte Carlo, as were the fractions of mixed events in each contribution. 

The parameters of the analytic functions describing the combinatorial background (N%, N\, 
M4(AM) and T^(t, Q) were fitted entirely from the data, with only the choice of functional forms mo- 
tivated by simulation. The shapes of the mass and proper time functions (including a small correlation 
between AM and t) are constrained by the same sign sample (which is almost entirely combinatorial 
background), and the opposite sign high AM region serves to normalise the number of combinato- 
rial background events in the low AM region. Since the combinatorial background is dominated by 
semileptonic b decays, it contains oscillating components in both opposite and same sign samples. 
The amplitude and frequency of this oscillation were constrained to be the same in both samples, but 
the fraction of unmixed events was allowed to be different, as supported by Monte Carlo simulation 
studies. 

The values of t-qo and Am^ were extracted by maximising the total likelihood given by equation ||. 
The values of t^o and Amd were allowed to vary, together with the number of signal events and 18 
auxiliary parameters describing the combinatorial background level, mass difference, time and mixing 
distributions. A result of 



r B o = 1.541 ±0.028 ps, 
Am d = 0.497 ± 0.024 ps" 1 
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Error Source 


A(r B o) (ps) 


A(Am d ) (ps" 1 ) 


B + lifetime 


0.004 


0.001 


B s lifetime 


0.002 


0.001 


Br(b -> B*+t~P) 


0.001 


0.001 


Br(b -> D*+d|* )- ) 


0.005 


0.002 


Br(b -> D*+X) 


0.001 


0.001 


b -> D**fl/X decays 


0.004 


0.009 


Other backgrounds 


0.006 


0.005 


b fragmentation 


0.005 


0.014 


D° decay multiplicity 


0.004 


0.003 


Flavour tagging offsets 


0.000 


0.001 


Flavour tagging mis-tag 


0.000 


0.009 


Tracking resolution 


0.015 


0.010 


Detector alignment 


0.003 


0.003 


Fit method 


0.012 


0.010 


Total 


0.023 


0.025 



Table 2: Summary of systematic errors on the measured values of t-qo and Am^. 



was obtained, where the errors are only statistical. The correlation between t-qo and Amd is —0.14. The 
fitted values of the other parameters were consistent with expectations from Monte Carlo simulations. 

The distributions of reconstructed proper time t for opposite and same sign data events with 
AM < 0.17 GeV, together with the fit results, are shown in Figure ||. The fractions of mixed events 
R, corrected for the dilution due to mis-tagged events, are shown for the same AM regions in Figure ||, 
where the oscillation of signal B° — * D* + £~v decays and the oscillation in the combinatorial background 
are clearly visible. The fit describes the data well, both in this mass region and in the high AM 
region dominated by combinatorial background. The discrepancy at high reconstructed proper times 
in the opposite sign sample (Figure |3](a)) is caused by imperfections in the resolution function (see 
Figure ||(f)) and is addressed systematic error. 



6 Systematic Errors 

Systematic errors arise from the uncertainties in the fit input parameters given in Table |l], the Monte 
Carlo modelling of the signal proper time resolution and backgrounds, the production flavour tagging 
performance and possible biases in the fitting method. The resulting systematic errors on the values 
of r B o and Am^ are summarised in Table || and described in more detail below. 

Input quantities: The various numerical fit inputs were each varied according to the errors given in 
Table [l] and the fit repeated to assess the resulting uncertainties. The B s oscillation frequency 
Am s was varied from the experimental lower limit of 9.1 ps -1 [22] to 50 ps -1 , with no significant 
effect on the fit results. 

b — > T>**£D decays: The systematic error due to the b — > ~D**£p background is dominated by un- 
certainties in the rate of such decays, determined from the measured Br(b — » D* + 7r~£pX) as 
discussed in Section |5[ To evaluate possible dependence in the b — ► D**£i) decay model, the 
default JETSET model used in the simulation was reweighted according to the calculations of 
Leibovich et al. [25] as used in |8], with negligible change in the fit results. The relative contribu- 



tions of B + and B decays to this sample are expected to be predicted well by isospin symmetry, 
but the B s contribution was scaled by 0.75 ±0.25 relative to the isospin prediction to account for 
possible SU(3) flavour violation effects, as in J|. Again, the resulting changes in the fit results 
were negligible. 
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Figure 3: Distributions of reconstructed proper time t for (a) opposite sign and (b) same sign events 
with AM < 0.17 GeV. The data are shown by the points with error bars and the expectation from the 
fit result by the histograms. The contributions from signal B° — > D* + £~is, resonant and combinatorial 
backgrounds are indicated. 
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reconstructed proper time (ps) 




reconstructed proper time (ps) 



Figure 4: Corrected mixed event fractions R vs. reconstructed proper time t for (a) opposite sign and 
(b) same sign events with AM < 0.17 GeV. The data are shown by the points with error bars and the 
expectation from the fit result by the solid line. 
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Other backgrounds: The backgrounds in fit source 3 from opposite sign b — » D* + tPX and b — > 
D* + Dg"X, and same sign D* + — > D°7r + , D° — > £ + X decays have characteristic lifetimes of 
1.87 ps and 1.42 ps respectively in the simulation. These values depend both on the true b 
hadron lifetimes and the mis-reconstruction of the apparent b hadron energy in these events, and 
were conservatively varied by ±0.1 ps to assess systematic errors. These backgrounds also have 
oscillating components due to the contributions from B° decays (typically 60% of the events). 
The systematic error was assessed by varying the B° fraction up to 100 % of the contribution 
from b hadrons, and repeating the fit. 

b fragmentation: The effect of uncertainties in the average b hadron energy (xe) = -^b/^beam was 
assessed in Monte Carlo simulation by reweighting the events so as to vary (xe) in the range 
0.702 ± 0.008 as recommended by the LEP electroweak working group [^q[ , and repeating the 
fit. This value is also consistent with a recent determination of {xe) = 0.714 ± 0.009 from 
SLD [27| |. The variation was implemented using the fragmentation functions of Peterson, Collins 
and Spiller, Kartvelishvili and the Lund group |p£fl , and the largest observed variations were 
taken as systematic errors. 

D° decay multiplicities: The b hadron energy reconstruction is sensitive to the B° decay multiplic- 
ity, which depends only on the D° decay for the B° decay channels of interest. The systematic 
error was assessed by varying separately the D° charged and 7r° decay multiplicities in the Monte 
Carlo simulation according to the measurements of Mark III [^9|. The branching ratio D° — ► K° 
or K° was also varied according to its uncertainty The resulting uncertainties on r B o and 
Anid from each variation were added in quadrature to determine the total systematic errors. 

Production flavour tagging: The jet and vertex charge distributions Q opp , Qs&me and Q v tx are not 
charge symmetric because of detector effects causing a difference in the rate and reconstruction 
of positive and negative tracks @, These offsets were measured directly from the data as 
in 0], and their statistical uncertainties contribute to the systematic error on Am^. 

The production flavour tagging mis-tag estimates must also be correct — i.e. that the values of 
Q correctly represent the mis-tag probabilities in the data. Excluding the events tagged by jet 
charge alone, the opposite flavour tags were checked in |Q] to a precision of ±2.6%, using events 
tagged in both hemispheres. The same technique has been used to check the jet charge Q op p, 
which was found to be correct to a precision of ±10%. The same hemisphere jet charge tag 
Qsame was also checked in Q to ±10%. These uncertainties were translated into systematic 
errors on Amd by scaling the values of Qt and Qm and repeating the fit. 

Tracking resolution: Uncertainties in the tracking detector resolution affect the time resolution for 
B° — > D* + £~v and B — > Y)**£~v events. The associated error was assessed in the simulation 
by applying a global 10 % degradation to the angular and impact parameter resolutions of all 



tracks, independently in the r-cj) and r-z planes as in [13], and repeating the fit. 



Detector alignment: The results are sensitive to the effective radial positions of the silicon detector 
wafers (both the positions of the detectors themselves and the positions of the charge collection 
regions within them), which are known to a precision ±20 fim from studies of cosmic ray events 
[13]. The resulting uncertainty was assessed in Monte Carlo simulation by displacing one or 



both microvertex detector barrels by 20 /mi and repeating the fit. 

Fit method: The entire fitting procedure was tested on a fully simulated Monte Carlo sample six 
times bigger than the data, with true values of r B o = 1.6 ps and Am^ = 0.437 ps -1 . The fit gave 
the results t-qq = 1.604 ± 0.012 ps and Amd = 0.440 ± 0.010 ps -1 , in both cases consistent with 
the true values. The statistical errors on the fitted values were taken as systematic errors to 
account for possible biases in the fit. Additionally, the large Monte Carlo sample was reweighted 
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to change the values of t b o and Anid, and the fit correctly recovered the modified values. To 
verify the correctness of the statistical errors returned by the fit, the large Monte Carlo sample 
was split into many subsamples, and the distribution of fitted results studied. The signal mass 
distribution Mi (AM) was varied to cover possible discrepancies in the Monte Carlo modelling 
of the signal mass spectrum (see Figure |l|). The data fit was repeated changing the minimum 
reconstructed time cut to — 1.5ps and the maximum to 10 ps (see Figure |3]( a)), and varying the 
maximum mass difference AM between 0.20 and 0.25 GeV. Finally, the data was divided up into 
5 subsamples according to the year of data taking. In all cases, consistent results were found, 
and no additional systematic error was assigned. 



7 Conclusions 

The lifetime and oscillation frequency of the B° meson have been measured using B° 
selected with an inclusive technique. The results are: 

r B o = 1.541 ±0.028 ±0.023 ps, 
Am d = 0.497 ± 0.024 ± 0.025 ps -1 



D* + £ v decays 




33, 34], with the world averages 



[ [22j| , and are in each case the most precise 



where in each case the first error is statistical and the second systematic. These results are consistent 
with other determinations of t b o [|]-||,[3l], 32 and Amj 
of r B o = 1.54 ± 0.03 ps and Am d = 0.470 ± 0.019 ps" 1 
determinations to date. 

The results presented here have been combined with the previous OPAL measurements of r B o 
using neutral secondary vertices Q and exclusively reconstructed B° — > D^iD decays |3l| and with 
measurements of Amd from single lepton (30] , dilepton |TyJ and exclusively reconstructed D* + decays 



dilepton [17] and exclusively reconstructed D* 
, taking into account statistical and systematic correlations. The results 

t b o = 1.538 ± 0.025 ± 0.023 ps, 
Am d = 0.479 ± 0.018 ± 0.015 ps -1 



were obtained, where again the first errors are statistical and the second systematic. 
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